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Abstract 

The BaR4R collaboration has nearly completed a program of precise measurements of the cross sections for the 
dominant channels of e^e^ — » hadrons from threshold to an energy of 3 - 5 GeV using the initial-state radiation (ISR) 
method, i.e. the measurement of the cross sections e^e" — > y hadrons with the energetic y detected at large angle 
to the beams. These data are used as input to vacuum polarization dispersion integrals, in particular the hadronic 
contribution to the muon g-2 anomaly. In addition to the recently measured tt^tt" cross section, giving the dominant 
contibution, many multihadronic channels have been investigated, with some recent examples presented here. We 
give preliminary results for the process e^e" — > K^K^{y) using 232 fb"' of data collected with the BABm detector at 
e^e" center-of-mass energies near 10.6 GeV. The lowest-order contribution to the hadronic vacuum polarization term 

K.K. LO 1 n 

in the muon magnetic anomaly is obtained for this channel: a^ ' - (22.95 ± 0.14(stat) ± 0.22(syst)) x 10" , which is 
about a factor of three more precise than the previous world average value. 
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1. Hadronic vacuum polarization and muon g-2 

An important part of the Standard Model prediction 
for the muon magnetic anomaly = (g - 2)/2, where 
g is the gyromagnetic ratio equal to 2 at lowest QED 
order, is given by hadronic vacuum polarization (HVP). 
In fact the dominant uncertainty in the prediction comes 
from the HVP contribution which is computed through 
a dispersion relation using the experimental information 
on the cross section for e^e" — > hadrons, as the relevant 
energy scale is too low for applying perturbative QCD. 
The HVP component is given by: 



had 



1 r°° 



(1) 



where K(s) is a QED kernel and ChajjonsC*) the bare 
cross section including final state radiation (FSR). 

2. The ISR method at BAK4R 

Unlike previous measurements that were done by en- 
ergy scans, the analyses presented here use the ISR 



method |IT|. The e^e^ — > X(y) cross section at the re- 
duced energy Vi' is deduced from the measured spec- 
trum of e^e^ — > X(y)yisR events produced at the center- 
of-mass (cm.) energy yfs. The reduced energy is re- 
lated to the energy £* of the ISR photon in the e^e^ 
cm. frame by i' = i(l - 2E'^/ ^/s). yfs' is equal to the 
mass of the final state X, including FSR photons. The 
ISR method follows from the relation 



dNx(y)yisR _ '^^I'sr 
d -/s' d yfs' 
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(2) 



where dL'^^^/d"^ is the effective ISR luminosity, exy 
is the full acceptance for the event sample, and cr^j^) is 
the 'bare' cross section for the process e^e" — > X(y) 
(including final-state radiative effects, but with lep- 
tonic and hadronic vacuum polarization contributions 
excluded). For precision measurements with the ISR 
method, the effective ISR luminosity is not taken from 
the theoretical radiator function, which describes the 
probability to emit an ISR photon of energy in a 
given angular acceptance, and the knowledge of the 
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e^e^ luminosity. Instead, it is determined from the mea- 
surement of the e^e" — » A(^//"(y)yisR spectrum on the 
same data sample, through a relation similar to Eq. Q 
where the e^e" — > ju^//" cross section is given by 
QED. In this way several systematic uncertainties can- 
cel. In particular, the measurement is mostly insensitive 
to higher order ISR corrections and other uncertainties 
affecting the hadron and muon channels equally. 

In the BaRar analyses the ISR photon is detected at 
large angle with £* > 3 GeV. This defines a topology 
where the ISR photon is back-to-back to the produced 
hadrons, thus providing high acceptance and better par- 
ticle identification (PID). Kinematic fits are used to re- 
ject backgrounds and improve mass resolution. A con- 
tinuous cross section measurement from threshold up to 
3-5 GeV is achieved, the upper range value depending 
on the background for each exclusive process. 



3. Overview of BaBar ISR results 

So far BaBar has published cross section results 
on 21 exclusive channels: tt^tt" ^ |3], 2{n*n') ||4l, 
;r+7r-;r", K^K-n+n-, K^R-ln^, 2(K^K-), KsK^'n^, 

K^K-jT^n-ri, 3(7r+;r"), 2(;r+7r-7r"), 2(7T+7r-)K^K-, (prj, 
0/"(98O), p-p, AA, AI^ and c.c, ijOr" |i51. 

Some final states are still under study: 7i^n^2jf, 
KsKl, KsKltt^tt-, KsK^tt'^tt'^, KsK^it^ri. Preliminary 
results on K^K are given in the next section. 

Some recently published results can be highlighted in 
view of their importance for g-2. The measurement of 
cr(e+e" — > ji^n^iy) is the most precise and complete 
for this process |2l |3]. The BaR4R ISR procedure is 
checked with the QED reaction cr(e+e" — > (y)). 
Both results are given in Fig. [T] showing agreement 
within 1.1% (dominated by the 0.94% uncertainty on 
the e^e^ luminosity) of the jUju(y) data with the next- 
to-leading-order (NLO) QED calculation. The jin{y) 
cross section measurement, being obtained from the ra- 
tio Tmlnii, does not rely on the e^e^ luminosity, thus 
providing a systematic uncertainty of only 0.5% in the 
dominant p region. 

Another important process is e^e^ — » 2{n^ji^) on 
which an analysis based on the full BaBm. integrated lu- 
minosity of 454 fb"' has been recently published Q. 
The cross section results are given in Fig. |2] together 
with those from previous experiments. Agreement is 
reasonable, the BaBm. results being both more precise 
and spanning the whole energy range of interest for 
HVP dispersion integrals. 
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Figure 1: (a) The ratio of the measured cross section for 
e*e^ — > jj^fi^yiy) to the NLO QED prediction. The band 
represents the best fit with the total uncertainty, (b) The mea- 
sured cross section for e*e^ — > n^K^iy) from 0.3 to 3 GeV. (c) 
Enlarged view of the p region in energy intervals of 2 MeV. 
The plotted errors are from the sum of the diagonal elements 
of the statistical and systematic covariance matrices from the 
unfolding procedure. 



4. Preliminary results on e'*'e -* K^K (7) 

4.1. Event selection and analysis 

The analysis is based on 232 fb"' of data collected 
by the BaBar detector |6j and follows closely the nn 
analysis |3|. Two-track ISR events are selected by re- 
quiring a photon with an energy E* > 3 GeV in the e^e^ 
cm. and polar angle with respect to the beam in the 
range [0.35-2.4] rad, and exactly two tracks of opposite 
charge, each with momentum p > I GeV/c and iden- 
tified as kaons with the DIRC Cerenkov detector and 
dE/dx. A K-ID efficiency of 80% is achieved, with fake 
rates at most 10% at the highest momentum. Kaon ID, 
as well as pion and muon mis-ID to kaons, are studied 
with high purity data samples. 

As for the analysis of the fi/jry and Tiny processes, the 
event definition is enlarged to include the radiation of 
one photon in addition to the already required ISR pho- 
ton. Two kinematic fits to e^e^ — > K^K^iy)yi5R are 
performed: (1) if an additional photon is detected in the 
EMC, with energy Ey > 20 MeV, it is used in a three- 
constraint (3C) fit, called 'FSR' fit (however the extra 
photon can be either from FSR or additional ISR); (2) 
a fit with the additional photon assumed to be emitted 
along the e'^ beam directions, called 2C ISR fit. This 
procedure allows the reconstruction of events at NLO, 
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Figure 2: The new BaBAR results (filled circles) on the e+e" — > 
2(ff^7r") cross section in comparison to previous experiments. 
Only statistical uncertainties are shown. 



therefore reducing the uncertainty resulting from the ne- 
glect of events with higher-order radiation. Each event 
is characterized by the two;^'p5j^ and;(f^gj^ values (except 
for the 12.5% with no extra measured photons), and the 
KK mass is obtained from the ISR fit if Xi^^ < AfpsR' 
and from the FSR fit in the reverse case. For the cross 
section measurement, the KK(j)yisR candidates are re- 
quired to satisfy InO^^j^j^ + 1) < 3. 

The overall acceptance is determined with a full sim- 
ulation, with corrections applied to account for observed 
data to MC differences. Through specific studies, the 
ratios of efficiencies are obtained in data and simula- 
tion for trigger, tracking, PID and selection, and ap- 
plied as mass-dependent corrections to the KK mass 
spectrum measured in data. Small corrections to the 
geometrical acceptance due to the assumed collinear 
additional ISR in the generator are obtained from the 
PHOKHARA code Q with fast simulation. 

Backgrounds stem primarily from other ISR events 
with extra ;r"'s or misidentified pions, muons,or pro- 
tons. Non-ISR qq represents the other important source 
of background, with an energetic photon from decay 
misidentified as the ISR photon. Estimation of back- 
grounds rely on PID measurements, rescaled simulation 
results using comparison with data in special samples, 
and the shape of the xjsK distribution. Background lev- 
els are negligible on the resonance, but increase at 
larger masses (about 2% at 1.1 GeV, 20% near 1.2 GeV 
fromp ^ 7T7T, 10% at 2 GeV, 50% at 5 GeV ). After spe- 
cific studies of the mass calibration and resolution, un- 
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Figure 3: The preliminary BaBAR e" — > K*K^{y) bare cross 
section including FSR. Systematic and statistical uncertainties 
are shown as diagonal elements of the total covariance matrix. 
The contributions of the Jj\p and 1^(25) resonances have been 
removed on this plot. 



folding of the background-subtracted mass spectrum is 
performed. Finally the effective luminosity is obtained 
with the /i/i(')') sample. 

4.2. Results 

The cr^x+x-(y)( V**) bare cross section including FSR 
is shown in Fig. [3] from threshold up to 5 GeV. It spans 
over a very large dynamical range (more than six or- 
ders of magnitude), and is dominated by the resonance 
close to threshold. Other structures are clearly visible at 
higher masses. 

The cross section in the (p region is given in Fig. |4j 
It is in fair agreement with previous CMD-2 IS] and 
SND f9l s results: the (p peak BABar cross section is 
about 5% (10%) higher than CMD-2 (SND), where the 
quoted systematic uncertainties are 0.7%, 2.2%, 7.1% 
for BaBar CMD-2, SND, respectively. There is also a 
small shift in mass (92 keV with CMD-2 and 65 keV 
with SND) consistent with the quoted mass calibration 
uncertainties, 110 keV (a preliminary conservative esti- 
mate) for BaBar and 80 keV for CMD-2). 

A fit of the charged kaon form factor measured by 
BaBar is performed taking into account the p and oj 
tails, as well as contributions from higher mass vector 
bosons. The following parameters are obtained: - 
(lO19.51±O.O2exp + O.llcai)MeV,r0 = (4.29+0.04exp + 
0.07re,soi) MeV, rfeB^.^- = (0.6344 ± 0.0059exp + 
0.0028fit + 0.0015cai) keV, where the last two uncertain- 
ties are from the form factor fit and the mass calibration. 
The precision of the last result is improved by a factor 
of two compared to the best value from CMD-2 jSJ. 
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Figure 4: The preliminary BAB/H e^e" — > K*K^{y) bare cross 
section including FSR in the (fi region. Data points from previ- 
ous CMD-2 and SND experiments are shown for comparison. 



5. The impact of BaBmi data on the g - 2 prediction 

The dominant nn HVP contribution to has received 
a lot of attention over the last 10 years, with discrepan- 
cies between experiments being partially resolved with 
time. Fig. |5] shows the results from all experiments as 
well as the determinations using t data corrected for 
isospin-breaking |10|. To keep results as independent as 
possible the comparison uses only data from the consid- 
ered experiment (complemented by world-average data 
in energy ranges not covered). It is seen that the ii4R4R 
result is the most precise (with CMD-2) and helps re- 
duce the tension between ee and r data. 

The precision of the new data presented here on 
the K^K^ and 2(ji^n^) channels also allow further 
progress. The contributions up to 1.8 GeV are 
^K-K-,hO ^ (22.95 ±0.14stat±0.22syst) x lO""' and 



a: 



= (13.64 ± 0.03s,at + 0.36sys,) X 10- 



10 



spectively to be compared with previous determi- 
nations im, (21.63 ±0.27smt + 0.68syst) x lO""' and 

(13.35 ± O.lOstat + 0.52syst) X lO""'. For other multi- 
hadronic channels the BfiB/H ISR results are by far the 
most accurate and complete. In addition the dynamics 
of each final state has been studied and found to be dom- 
inated by resonances which have been identified. This 
is important because it is possible in this way to derive 
some cross section estimates llT2l for final states which 
are difficult to measure (> lif*) using the known branch- 
ing fractions of these resonances. This detailed infor- 
mation has been used in the most recent estimate of a^, 
showing a 3.6cr discrepancy fV2\ with the direct mea- 
surement 1 13 1 . Another estimate L14j uses also ii4R4R 
cross section data. 
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Figure 5: The nn HVP contribution to obtained from t de- 
cays with isospin-breaking corrections (top) and e^e" data 
(bottom) fTOlflTl . 
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